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Resume 

We show how to cope with the acoustic identification of poroelastic materials when the specimen is in 
the form of a cylinder We apply our formulation, based on the Biot model, approximated by the equivalent 
elastic solid model, to a long bone-like or borehole sample specimen probed by low frequency sound. 

1 Introduction 

It has become fairly common fP[, fT\, f4], f5], f6], fT\, fsl, f9], fTP[, fTT\, fl3], (14], fB), 
to acoustically identify the structural/material properties, and/or the phase velocity/attenuation in poroelas- 
tic (e.g., biological or geophysical) materials by processing data relative to the reflected and/or transmitted 
pulses of a slab-like specimen of the material in response to an incident plane- wave pulse. In most of these 
studies, the model of the medium is that of a fluid or (less often) viscoelastic solid. It is not often that mate- 
rials, especially those of biological 1 16 1 or geophysical 1 17 1 nature, present themselves in the form of slabs, 
plates or layers, nor is it judicious to cut them (an operation that is rarely accurate and which often modifies 
the physical properties) to fit this shape. We show herein how to cope with the acoustic identification of 
these materials, in the form of cylindrical specimens |fl61, fll^, whether they are considered to be fluid-like 
or elastic solid-like. Work on the fully-poroelastic (in the sense of Biot) case is in progress. 

More specifically, this investigation is concerned with the reconstruction of the material constants = 
(k^ the complex compressibility), fj,^ (shear modulus), (density) of an elastic solid-like (assumed to be 
of this nature even if the target is poroelastic in the sense of Biot), almost-circular cylinder, modeled as a 
circular cylinder of radius a close to the average radius of the almost-circular cylindrical target. The latter 
is immersed in a fluid-like host wherein propagates plane-wave like acoustic probe radiation. The action of 
this wave on the target results in a scattered acoustic wavefield which serves to reveal the material properties 
of the target. 

The wavelength (A^ — 27r/fc° ) of the probe radiation is assumed to be much larger (x^ '■— k^a << 1) than 
a. A perturbation analysis is shown to enable an explicit reconstruction of and of a linear combination of 

This reconstruction technique relies on the a priori knowledge of a ; it is assumed herein that this parameter 
is known. 



2 Physical configuration 

The scattering body is an infinite cylinder whose generators are parallel to the z axis in the cylindrical 
coordinate system (r, 6, z). The intersection of the cylinder, within which is located the origin O, with the 
xOy plane defines : 

i) the boundary curve T = {r — f{9);0 < 9 < 27r}, with / a continuous, single-valued function of 9; 
further on, we shall take F to be a circle, i.e., f{9) = a, with a its radius, close to the average value of ri{9), 

ii) the bounded (inner) region (i.e., the one occupied by the body in its cross-section plane) ili = {r < 

r]{9); 0<9 < 27r}, 

iii) the unbounded (outer) region flo = {r > ri{9); Q < 9 < 2tt}. 

It is assumed that VIq is filled with a linear, homogeneous, inviscid fluid and ili by a linear, macroscopically- 
homogeneous, isotropic, porous medium which will subsequently be associated with a linear, homoge- 
neous, isotropic, time-invariant elastic solid medium. 

The material constants of are assumed to be known. Those of Af ^ are unknown and are to be recovered 
by the technique described hereafter The latter relies on probing the cylinder (from the outside) by a plane 
acoustic wave whose wavevector lies in the xOy plane. 



3 Ingredients of the Biot theory 

We give the ingredients of the basic Biot theory of biphasic (soHd/fluid) porous media 1181 . II19I . 



3.1 Conservation of momentum relations 

In the absence of appHed body forces, the conservation of momentum relations take the form 

r/F 

V • cr - pu - PfVf,tt = , Vp + PfUtt + mw^tt H = , 

K 



(1) 



wherein : 



- i t designates a first-order partial derivative with respect to time t, and f tt := {f,t),t, 

- (T is the total stress tensor in the porous fluid-saturated medium, 

- ps the density of the solid component 

- Pf the density of the viscous fluid filling the (interconnected) pores, 

- p the bulk density of the porous medium, such that p — {1 — (p)ps + (ppf, with (p the porosity (volume 
fraction of fluid relative to total volume in a representative volume element), 

- a the stress tensor, 

- u the displacement vector of the solid particle component, 

- U the fluid particle displacement vector, 

- w the relative displacement vector of the fluid particle relative to the sohd particle (both particles assumed 
to occupy the same point) defined as w := (/)(U — u), 

- p the pressure in the fluid component of the porous medium, 

- rj the viscosity of this fluid, 

- K the (low-frequency) permeability, 

- m — the virtual mass, 

(p 

- a the tortuosity (in l20l a is termed the virtual mass coefficient or structure factor), which, in 1211 . is 
related to the porosity by a = 1 + r ("'^^ (wherein r is a constant with a value between and 1), 

- F{t) is a linear integral convolution operator with respect to time which, in the frequency domain, becomes 
a frequency-dependent multiplier F{uj), implying frequency-dependent permeability, i.e., k{u!) = -p^, 

wherein k{lu) is the so-called dynamic permeability, and F is designed so that lim^^o F{u!) = 1. 



3.2 Constitutive relations 

Biot's constitutive relations f 181, 11221 linearly relate the total stress and fluid pressure to the (isotropic) 
solid and fluid displacement spatial derivatives via 



a ^2pe+[{H -2p)e-CC]I , p = Af ( - Ce 
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wherein : 



e is the strain tensor £ = 5 



^ fVu + Vu^) 

I the unit tensor 

- p the shear modulus (rigidity) of the saturated solid, 

- K = {X + Kg the bulk moduli of the dry (i.e., drained) solid matrix and solid grain material respec- 
tively (note that Ac= ^1 — -^^C + Ais the Lame constant of the saturated solid), 

- Kf the fluid bulk modulus. 



3.3 Equations of motion in terms of u and w 



We assume henceforth that all the material parameters are constants with respect to position, i.e., the 
medium is macroscopically homogeneous. The conservation of momentum equations and constitutive rela- 
tions are employed in such a way as to eliminate the pressure and stress tensor so as to obtain : 

/iV2u+(i?-^)Ve-CVC = pU(t+p/w,tt , CVe-MVC = P/u„+mw,tt + -Fw,t , (5) 

K 

which is the vectorial form of the Biot wave equations as given by Yamamoto 1 19] for a macroscopically- 
homogeneous porous medium. 

4 Choice of an approximate model to describe wave propagation in a 
porous medium 

The fundamental difficulty with the Biot theory is twofold : i) two coupled (vectorial) wave equations 
have to be solved simultaneously, and ii) many material parameters have to be recovered in the inverse 
problem context. This is why the traditional approach (notably in the underwater acoustics community) has 
been to reduce this model to a simpler one (with fewer material parameters) in which only one (vectorial or 
scalar) wave equation has to be solved. 

4.1 Equivalent elastic solid model (EESM) obtained from the limit 0^0 

The introduction of the second equation of (|5} into the first gives 

AiV^u +(h - ii--\ve+ (-C + — ) VC + (-P + ^) Utt + —F^v.t = , (6) 
\ m J \ m J \ m / mn 

Recalling the definitions of m and w, we conclude that lim^^o = and lim^^o w = 0, so that in 

the limit ^ 0, (|6} becomes 

^V^u + (Ac + VV • u - putt = , (7) 

which is simply the Navier wave equation in a non-dissipative, linear, homogeneous, isotropic solid in which 
the material parameters are Ac, /i, p — Ps- 

This equation forms the basis of the so-called equivalent elastic solid model for wave propagation in poroe- 
lastic media, often employed for the evaluation of transmission loss of sound over sediment layers on sea 
bottoms L23J. 124J. 

4.2 Equivalent viscoelastic solid model (EVSM) for wave propagation in poroelastic 
media 

We now consider another approximation of the Biot wave equations which leads to what has been termed 
the equivalent viscoelastic solid model of propagation in poroelastic media l,25ilL22J . 



The basic idea is to reduce the Biot model to only one wave equation, while retaining the loss mecha- 
nism inherent in this model (contrary to what is done in the equivalent elastic solid model in which the loss 
mechanism is abolished in the limit — * 0). 

Let us return to the definition of w, i.e., w = (/)(U — u) in which we assume 0U << (f)U. Thus, ne- 
glecting the terms in (/)U in (|6}, we get 



H - p C0 H - 
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VV-u+(-p+^)u„-^Fut=0. (8) 

TO / rriK 



This wave equation for the displacement in the solid component is similar to the Navier equation, with the 
exception that the wavenumber is now complex, its real part not being equal to that of the Navier equation 
wavenumber, and its imaginary part being conditioned by 



4.3 Equivalent fluid model (EFM) 

The equivalent fluid model (EFM) is appropriate when the fluid is light (e.g., a gas such as air) and the 
solid skeleton is therefore relatively immobile (i.e., rigid). This model has been employed, even when the 
fluid is not light, notably in the underwater acoustics community 1261 . 1231 . 1191 . 

There exist various versions of the EFM, several of which are described and compared in DepoUier et 
al. 1271 . The one we shall consider herein is a simplified version of the model offered in II II . 

No restrictions are introduced concerning the porosity, but, for the sake of simplicity, the fluid viscosity 
is assumed to be nil. The fundamental assumption is that the solid component is rigid, i.e. 

X, = c5o, , u = -^u,tt = 0. (9) 

If, in addition, the medium is macroscopically-homogeneous and time-invariant (i.e., pf, a and Kf are 
constants with respect to position and t) then the Biot system of equations reduce to 

-V^p-p/aV-Utt = , p^u + Kf\J,u^O, (10) 
a linear combination of which yields 

V'P-^P.« = 0, (11) 

which is the wave equation in the equivalent fluid. 
Note that a — 1 for a homogeneous fluid. 

4.4 Our choice of approximate model 

Since all three approximate models reduce to the same type of wave equation (either vectorial or scalar), 
a generic choice-( vectorial in nature) adopted herein-is the EESM. 

This means that we replace, by thought, the porous medium cyUnder immersed in an inviscid fluid by an 
elastic solid cylinder immersed in the same fluid. 

5 Mathematical description of the problem 

5.1 Preliminaries 

Due to the invariance of the cylinder and incident wavefield with respect to z, the incident and scat- 
tered fields are also invariant with respect to z. 

Let designate pressure in ilo ; due to this invariance, p^ = p(a;, y, t) — p^{x, y, t) + p'^[x, y, t), wherein 
p'^ is the diffracted pressure in f^o. 

For the same reason, the total displacement wavefield u in 57i is of the form u = u(a;, y, t). 

The analysis is carried out in the space-frequency domain via p°'*'''(x, t) — p^'^''^{x, uj) exp{—iLijt)duj 
and (x, t) ~ (x, uj) exp{—iujt)duj wherein x — {x, y). 

Henceforth, it is implicit that * (meaning p^, or p*^) and are functions of (x, lo). 

5.2 Governing equations 

The incident plane wave is 

p' ^ S{uj)exp{-ikorcoB{e -e')) (12) 
(wherein S{uj) is the amplitude spectrum, 6'* the incident angle in the xOy plane, and kF' — iojtP, with 



pO,t,d satisfy the frequency-domain pressure wave (Helmholtz) equation 

(V^ + (fc°)2)/-^^^ = in 1^0, (13) 

and the radiation condition 

- i^V = o(r"i/2) : r ^ oo , e [0, 27r[ . (14) 

satisfies the frequency-domain elastic wave equation 

/uV^u^ + (A^ +/ii)V • Vu^ +a;2ui = in f^i , (15) 

and the boundedness condition 

<oo in ill . (16) 

Let v designate the unit outward-pointing (from il^) unit normal vector, — cr^ ■ v the traction. Then the 
transmission boundary conditions are : 

T° - = , u° • 1/ - ■ IV = on r . (17) 

5.3 Reduction of the elastic solid wave equation to two Helmholtz equations 

The use of the Helmholtz decomposition 

= V(^i + V X -0^ (18) 
enables M5\ to be reduced to the two (one scalar, the other vectorial) Helmholtz equations 

(v2 + (fci)2)^i^0 , (V2 + (4)2)^^1=0, (19) 

wherein 

Recalling that the fields p and u do not depend on z enables (I13> and (I19> to be cast into the cylindrical 
coordinate forms : 

+ r- V'. + r-''-,p]gg + {k^TfiP' = . (21) 

The gauge condition V • ■0^ =0 and the absence of shear stress in the fluid imply that i^l — i^}, — 0. In 
the cylindrical coordinate system, the traction and normal component of displacement continuity conditions 
reduce to : 



^^j^p",-</)i,-r-Vi,, = 0. (22) 



5.4 Field representations 

The incident pressure field satisfies the periodicity condition pUr, ~0 + 2(?*, lo) = p^{r, 9, uj) and the 
first of the Helmholtz equations in i2l\ so that (also on account of il2\ 

oo 

P" =Y1 «mem^m(fci?-) COS to(61 - 6') , with a,„ = S{uj)e-'"'^ , (23) 

m=0 

wherein J,„( ) is the m-th order Bessel function and eq ~ 1 , em>o = 2. 

The periodicity of entails p'^{r, —9 + 29'', uj) — p'^{r, 9, uj), so that on account of the first of the Helmholtz 
equations in i21i and the radiation condition 

oo 

/ = J2 brne„,H^{klr) cos m{9 - 9') , (24) 
wherein H„i{ ) ~ H^} ( ) is the m-th order Hankel function of the first kind. 

The periodicity of p' and p'^ also entails 0^ (r, —9 + 20*, lS) — cj)^ (r, 9, lo), so that on account of the second 
of the Helmholtz equations in \2\\ and the boundedness condition 

oo 

0^ = ^ c,ne,nJniiklr) COS m(6' - 9') . (25) 

m=0 

By means of any one of the transmission conditions J22> . and on account of the periodicity conditions 
satisfied by p\ p"^, and c/)^, it is found that -0] obeys the relation '0](r, —9 + 29\ uj) — —4>l{r, 9, lo) so that 
on account of the third of the Helmholtz equations in \2\l and the boundedness condition 

oo 

■0] = ^ dmf^mJmikxr) smm{9 - 9') . (26) 

m=0 

6 Use of the transmission boundary conditions to obtain {bm}, {cm}^ 

{dm} 

We employ the orthogonality relations 

COS m{9~ 9') COS n{9- 9')— , / smm{9 -9') sin n{9- 9')— = 5., ^ 



7 I _ ^ - ^ ^mn 7 

TT em Je' TI" ^ 

(27) 

wherein Smm = 1 , <5mn/m = to obtain ; 

Poqo = ro , (28) 

in which 

fxlHoixl) -vlxlMxl) \ /bo\ f-aoxlMxl: 

Po = , qo = , ro = 

\~Hoixl) uiJo(xi)-2M^(fci)2jo(xi)y VoJ \ aoMxD 

(29) 

and 

P«qn = Tn ; 1,2, ... , (30) 



in which 
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with = fc^a, Xt — k^o,, ''^i = •^"' (^1)^^ = dg/dc;, g — cPg/d<,'^, and for n = 1,2, 

In principle, the matrix equations (I28> and ( I30> enable to determine the unknown coefficient vectors qo 
and q„ ; n = 1, 2, and thus to solve the forward-scattering problem, notably for the prediction of the 
scattered pressure field pd in the host fluid. 

Rather than do this, and since we are more interested, in the present context, in solving the inverse-scattering 
problem of the reconstruction of A^, /z^, p^, we adopt a different strategy for determining qo and q„ ; n = 
1,2,.... 



7 Low-frequency approximation of the solution of the forward-scattering 
problem in the region outside of the body 



We now define this new strategy. 
We first note that : 

xl^T^xl = -tx , XT^^xl = -tx, (33) 

wherein % :— x^. We assume that x is small enough (i.e., < x << 1) for it to be true that < xi << 1 
and < Xt < < 1' ^'^'^ employ a perturbation scheme, based on the smallness of x (which, to the very least, 
implies very low frequencies uj and/or small cylinder radius a), to solve the matrix equations. 

Thus, the arguments of all the Bessel and Hankel functions appearing in the expression of P„ are small, 
which fact authorizes use to be made of the small-argument asymptotic forms 

Jr.{0 - (0 " ; - I me , , H^iX) ^ -'-^^ (f ) ; c - , = 0, 1, ... , 

(34) 

To do this in a systematic manner, we expand P™, q^, and in series of powers of x : 

j=0 n=0 



l'm(x) — ^ ] I'm X i I'm J\'Q^^"iix) \ ^—q i (36) 

which, after introduction into the matrix equation PmQm — I'm, yields (after comparison of powers of 
X) X]L=o Pm~"^Qm'' = I'm ! ^ = 0, 1, 2, wWch defines the recursive scheme for the determination 



of qL : 



n(0)^fp(0)A ^„(0) / ^ (-p{0)\ ^ 

Mm \ m I m i Vim \ m I 



l-l 



r(0 _ V p('-")n(") 



1 = 1,2,. 



After a series of algebraic manipulations, the following asymptotic form of &,„ is found : 



(37) 



(38) 



wherein 



^0 = ao ( — 



r,(2) _ 
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Ai+A^i 

so that the diffracted pressure field in the host fluid becomes (to second order in x) 



(39) 



(40) 



8 Explicit resolution of the inverse problem : recovery of the material 
parameters of the specimen 

By making use of ( I40> and the first of the orthogonality relations \21\ we find 



Ai bl 



(2) 



4i 



p'^{b,e,uj) COS m{e - 9')— ; m = 0, 1 , (41) 

TT 



which signifies that Bq and Bi can be obtained from integrals involving the measured diffracted pressure 
field data (for all angles 9) on a circle of radius b. Once these two coefficients are found, A^, /i^ and can, 
in principle, be obtained from i39\ . i.e.. 



Bo 



X + fi-1 
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1 + PJ 



wherein 



These relations : 
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(42) 



(43) 



- show that the field is a nonlinear function of the material parameters, and 

- apply equally-well to the fluid model of the specimen (i.e., p = 0). 

Moreover, since only two pieces of data (i.e., Bq and Bi, as expressed by \A\V ) are available in this (second) 
order of approximation (in x), only one (p), and a linear combination (A + p) of the other two of the three 
material parameters can be recovered via (I42t . 



More specifically, it is found that 



A + p = 



1 



P = 



1-Bi 
l + Bi' 



which underlines the fact that the proposed technique enables 

- an explicit (partial) solution of the inverse problem of the reconstruction of p and "^|^u 
when the material of the specimen is modeled by an elastic solid, and 
an explicit (complete) solution of the inverse probh 
the material of the specimen is modeled by a fluid. 



(44) 



P 



an explicit (complete) solution of the inverse problem of the reconstruction of p and 7^42 = - when 



9 A manner for obtaining the complete solution of the inverse prob- 
lem for the elastic solid model of the specimen 

The previous analysis showed that Bn = 0;n = 2,3,....to second order in x, which is the reason why 
only two pieces of data are available to reconstruct the three unknown material parameters p, A and /x. 

It can be shown that by carrying out the perturbation analysis to fourth order in x, not only Bq, Bi, but 
also B2 are non- vanishing, so that by this means one disposes of the three pieces of data necessary to recon- 
struct the three unknown parameters. 

The only difficulty with this procedure is that the relations between Bq, Bi, B2 and p, A, p are much 
more complicated than previously, which fact makes it impossible to obtain expUcit algebraic expressions 
for p, A, p in terms of Bq, B\, B2. 

It turns out that one is faced with the problem of solving a system of three nonhnear equations in three 
unknowns. The procedure for solving this system is advantageously initialized via the second-order-in-x 
approximations of\ + p and p. 

10 Conclusion 

We have shown that the low-frequency perturbation scheme enables an explicit reconstruction of two 
(A^ + p and p) of the three material parameters p, \, pofa poroelastic cylindrical specimen, modeled as an 
equivalent elastic solid circular cylinder. 

Two single-frequency pieces of data, both involving the scattered pressure field on a complete circle around 
the cylinder, are required for this procedure. 

Data obtained at other frequencies can be employed to reconstruct the frequency-dependent complex pa- 
rameter A^(w) + p{uj) of cylindrical specimens modeled as viscoelastic materials. 

Above all, the same type of perturbation analysis can be employed to obtain appropriate combinations of 
the Biot parameters when the cylindrical specimen is modeled as a fuUy-poroelastic material (in the sense 
ofBiot). 
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